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Abstract 

We show how a boundary state different from the (1, 1) Cardy state may be realised in 
the (p,p + 1) minimal string by the introduction of an auxiliary matrix into the standard 
two hermitian matrix model. This boundary is a natural generalisation of the free spin 
boundary state in the Ising model. The resolvent for the auxiliary matrix is computed using 
an extension of the saddle-point method of Zinn- Justin to the case of non-identical potentials. 
The structure of the saddle-point equations result in a Seiberg-Shih like relation between the 
boundary states which is valid away from the continuum limit, in addition to an expression 
for the spectral curve of the free spin boundary state. We then show how the technique may 
be used to analyse boundary states corresponding to a boundary magnetic field, thereby 
allowing us to generalise the work of Carroll et al. on the boundary renormalisation flow of 
the Ising model, to any (p,p + 1) model. 
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1 Introduction 

Defining quantum gravity by a path integral over geometries has been a technically difficult 
idea to realise. One approach to making such an idea precise is known as Dynamical Trian- 
gulation (DT), which consists of approximating each geometry by a discrete geometry built 
from gluing together elementary building blocks, most often triangles, or higher dimensional 
analogues. The program of DT is most often carried out in two dimension, however it was 
hoped for a number of years that DT would provide a definition for quantum gravity in 
dimensions greater than two. This unfortunately did not turn out to be the case (for an 
overview consult [U |2]) and attention has since shifted to the use of a refined method known 
as CDT, introduced by Ambj0rn and Loll [3]. 

The theory of DT in two dimensions is nevertheless interesting, as the sum over two di- 
mensional random surfaces may be reinterpreted as the worldsheet theory for the Polyakov 
string [I]. The choice of matter coupled to gravity is then interpreted as describing the 
background in which the string is embedded. Furthermore, in two dimensions there ex- 
ist analytical methods for studying the theory. Firstly, the continuum path integral over 
geometry can be shown to be equivalent to quantum Liouville theory, which is integrable. 
Secondly, the sum over the discrete planar graphs required in the DT approach may be 
explicitly performed via matrix integrals, thereby allowing an analytic investigation of the 
properties of DT (for a review of these topics consult [2J EJ) • Thus we have two comple- 
mentary approaches in two dimensions. An important step in relating these two approaches 
was the work by Kazakov and Boulatov which provided a two-matrix model formulation of 
the Ising model coupled to DT [6J [7] (see also [H [9]) and the seminal works of Knizhnik, 
Polyakov and Zamolodchikov (KPZ) [10] and David [UJ and Distler and Kawai [12]. Addi- 
tional results, in particular, the formulation of the double scaling limit, came from the works 
of Brezin and Kazakov [13] , Douglas and Shenker [14] and Gross and Migdal |15j . 

In general, adding a statistical mechanical system to the DT triangulations or equivalently 
the addition of a minimal model (p,q), as introduced by Belavin, Polyakov and Zamolod- 
chikov [16] , to the Liouville theory leads to what is know as a (p, q) minimal string [T7 ] [T5 | [T9"] . 
This point of view was responsible for much of the early interest in DT, and led to the discov- 
ery of many connections to integrable hierarchies [2] . In this perspective boundaries play an 
important role as they correspond to branes. In the continuum formulation of Liouville the- 
ory coupled to a CFT, two forms of boundary states for the gravitational degrees of freedom 
have been identified; Fateev-Zamolodchikov-Zamolodchikov-Teschner (FZZT) [20, 21] and 
Zamolodchikov-Zamolodchikov (ZZ) [22] which depend on a real parameter a and integer 
valued parameters n and m respectively. When the FZZT boundary condition is tensored 
with a Cardy state for the minimal model degrees of freedom, one obtains a brane known as 
a FZZT (k,l) brane, denoted by |er; k, I) . On the DT side however this situation is not as 
clear. A two hermitian matrix model of the form, 



is required in order for the scaling limit to realise a (p, q) CFT interacting with Liouville 
theory. However, this matrix model has only two observables that correspond to boundary 




(1) 



2 



states; the resolvents Wx(x) = (-^Tr^-^r) and Wy(y) 



x 




These resolvents 



act as generating functions for boundaries of finite size. The existence of only two obvious 
boundary states is in contrast to the (p — l)(g — l)/2 boundary states one can obtain by 
tensoring a Cardy state with the FZZT boundary condition. 

The main suggestion for resolving this issue was given in [23J , in which it was noted that 
if the boundary state is inserted into a spherical worldsheet (thereby creating a disc), then 
the following identification of states holds, 



which is often referred to as the Seiberg-Shih relation. Since it was noted that the resolvent 
for X produced the disc amplitude corresponding to insertion of the state \a; 1, 1) it was 
therefore claimed there was no contradiction between the DT and Liouville approaches (see 
also [24J). This was later challenged by a number of authors [251 EE], who observed that 
such an identification of states no longer held when the state was inserted into a worldsheet 
of more complicated topology. However, the issue of whether this is indeed the case is yet 
to be entirely settled [27]. Despite this, we feel it is still interesting to understand if the 
other boundary states may be computed directly from the matrix model and secondly to 
understand where the Seiberg-Shih relation appears in the matrix model in the planar limit. 
These questions have been partially addressed by [28] in which they implicitly argue the 
Seiberg-Shih relations should hold away form the scaling limit. However, they do not address 
the construction of the spectral curve nor how to relate their construction to the topological 
recursion relations of Eynard and Orantin [29]. Finally, the physical interpretation of the 
boundary states constructed therein is less clear. 

In this paper we provide a different perspective on the solution to these questions, which 
generalises the work of [30, 26J in which a boundary magnetic field applied to an Ising model 
on a fluctuating lattice was studied. We show how one can use saddle point equations to 
study a boundary magnetic field applied to any (p,p + 1) minimal string. In performing 
this study we will see Seiberg-Shih like equations arise naturally at the discrete level before 
a scaling limit is taken. Such relations naturally generate algebraic curves to which the 
construction of [29 could be applied. For a non-zero boundary magnetic field we are able 
to extend the results of [30] [26] , in which fixed points of the boundary renormalisation flow 
were found, corresponding to fixed and free spin boundary states, to all (p,p + 1) minimal 
strings. 

2 Free spin boundary for the (p,p+ 1) minimal string 

Recall the two hermitian matrix model, 



fc-i 



1-1 




(2) 




(3) 
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where the potential V\ is of order p and V2 is of order q. The (p, q) minimal model coupled 
to Liouville gravity is the continuum theory corresponding to the scaling limit of this two 
hermitian matrix model about its highest order multi-critical point. Furthermore the (1, 1) 
FZZT brane may be realised in the matrix model by the resolvent for X. The resolvent for 
Y corresponds to the dual FZZT brane [23] • These resolvents may be computed using the 
master loop equation or spectral curve 



E(x, y) = {V[{x) - y) (V^y) - x) + gP(x, y) = 0, (4) 

where the resolvent Wx(x) = -(U'(x) — y(x)) and P(x,y) is a polynomial in x and y with 
undetermined coefficients that require extra analytic data (such as the one-cut assumption) 
to be determined. 

The (p,p + 1) minimal models coupled to Liouville gravity may be obtained in a second 
way, distinct from that given above, by having V\ and V2 identical and both of order p. 
This second way of realising the (p,p + 1) models will be of central importance in this 
paper. In particular for the case p = 3, which corresponds to the Ising model, there is a 
clear interpretation of the matrix degrees of freedom as Ising spins; X and Y vertices in a 
Feynman diagram correspond to up and down spins respectively. With this interpretation 
it is easy to construct a representation of the Cardy state corresponding to free spin [32] ; it 

corresponds to the X + Y resolvent, Wx+y{z) = (j^^ z -qc+y ) 

For p > 3 we expect the X + Y resolvent to again flow to a conformal boundary condition 
in the scaling limit that differs from the identity boundary condition obtained from the X 
resolvent. We will refer to this boundary condition as the free spin boundary state. 

The problem of calculating such a resolvent for the case p = 3 can be solved by simply 
changing variables in fl3]) to S = X + Y and A = X — Y. This converts ([3]) to a 0(1) model 
and the X + Y resolvent into a standard resolvent of the matrix S, which may then be 
computed for example via loop equations [32l [33] . However, this method does not generalise 
in an easy way to larger p. Instead we will now present an alternative calculation of the free 
spin boundary condition of the Ising model which is ripe for generalisation. 

First we consider Wx, this enforces a fixed spin boundary condition by ensuring that 
a spin residing in a face touching the boundary sees an up spin on the other side of the 
boundary edge as shown in Figure [H It does not however put any restriction on the spins 
appearing within the triangulation. In contrast, we now proceed by making the assumption 
that a boundary condition can equally well be introduced by imposing a constraint on the 
spins residing in faces touching the boundary as shown in Figure [2J This form of constraint 
can be imposed by introducing another matrix into the matrix model action which when 
integrated out yields the original matrix model. For example consider the following model£] 



Z = J [alM}[dX][dY] exp f-^Tr^pf) + V 2 (Y) + ^M 2 - XY - XM}^j 



(5) 



Note that M only has the bulk effect of renormalising the X propagator. However if we 
compute the resolvent of M it will produce Feynman diagrams, as show in Figure |2l in 



1 In this section we always denote an arbitrary potential by Vi without necessarily meaning they are the 
same potentials in each expression. 
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Figure 1: An example of a Feynman 
diagram contributing to (TrX k y In 
the bulk we have drawn the dual graph, 
with faces coloured dark and light grey 
corresponding to X and Y vertices re- 
spectively. The external X k vertex en- 
forces an up spin to reside on the ex- 
ternal side of each boundary link. 



Figure 2: An example of a Feynman 
diagram contributing to (TrM k ^ in the 
case M only couples to X. The M — > 
X propagator is denoted by a black-to- 
dark grey line. We see that the exter- 
nal M k vertex enforces an up spin to 
reside on each of the faces touching the 
boundary. 



which faces lying on the boundary may only contain up spins. In this example we expect 
the M resolvent will be equal to the X resolvent, up to non- universal terms. 

Generalising this idea, we see we can introduce an extra matrix M, which by coupling 
equally to X and Y, will have a resolvent computed by summing over all configurations of 
spins in the faces bordering the boundary as shown in Figure [3j Explicitly we have, 

Z = J [dM}[dX][dY] exp (-NTrlVtiX) + V 2 (Y) + ^M 2 - XY - (X + Y)M]^j , (6) 

where we have absorbed the g dependence in to A. We expect that resolvent for M will flow 
to the same boundary condition as the X + Y resolvent, in the scaling limit. 

The matrix model in ([6]) corresponds to a 3-state Potts model with non-equal potentials. 
The Potts model has been solved via loop equations but only for all potentials identical 
[31]. However, there exists an alternate solution of the Potts model based on saddle-point 
methods which can be generalised to allow for non-identical potentials. Firstly, note that 
the above model may be diagonalised by the trick of Daul [3S]; introduce a new matrix A 
such that the model may be written as, 

Z = j [dM][dA}[dX][dY}exp^-NTi[V 1 (X) + V 2 (Y) + ^M 2 -(X + Y + M)A + ^A 2 ]^j, (7) 

where V\(x) = Y^k=i ^t xk + Jg x2 an< ^ ^( x ) = J2k=i ~k xk + Jg x2 ' Note that g dependence 
is now absorbed in 7. We now may integrate out the angular degrees of freedom using the 
Itzykson-Zuber integral applied to each of the XA, Y A and MA terms. The resulting model 
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Figure 3: By allowing M to couple to X and Y equally, we have that (TiM k ^ will 
be computed by summing over diagrams in which up and down spins appear with 
equal weight on the boundary. Here the black-to-light-grey line represent M — > Y 
propagators. 

falls into a general class of models of the form, 

r Q / Q \ 

Z= [dA\ Y[[dB®} eX p -iVTr[^ V5(fl (<) ) + V A (A) - A • (8) 

In a work by Zinn- Justin [36] the above model was solved for Q < 5 in the case Vi = V 
for all i, via a saddle point method. For our purposes we need to solve this class of models 
with unequal potential, fortunately the generalisation to the case of differing Vi is simple. 
First we review the method of Zinn- Justin [35] for the case Vi = V for all i, for which the 
integral may be written as, 

Z = j [dA}e- NTrVA ^ {E[A]f . (9) 

where 

E[A]= f[dB}e NTr ^ v ^ +AB l (10) 



If we compute the saddle point equations for the eigenvalues tij of A, we will encounter the 
function ^- log H [A] . Following [36] we introduce the resolvent for A, Wa(cl) and a new 
function b(a) defined such that for large iV we may write, 

1 (> logE[A] = t(a J )-W A (a J ), (11) 



Nddj 



where f(z) = ~(f(z+i0)+f(z—i0)). Note that the potential V(B) plays no role in this saddle 
point equation, the resolvent Wa arises from the variation of the Vandemont determinant 
and b(a) is due to the variation of the AB interaction term in the action. Varying HL4] now 
with respect to the eigenvalues of B, b^, one obtains, 

^ l0gH[A] = m + ^ B{hj) ~ V ' {h ^ (12) 
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where we have introduced a new function a (6) to account for the variation of the AB term 
with respect to the eigenvalues of B. The crucial result of [36] is that 6(a) possess the 
following properties, 

• The function 6(a) has the same cut as the physical cut of the Wa resolvent, i.e. the 
location of the cut is the same and the difference across the cut is identical for the two 
functions. 

• The function b(a) is the functional inverse of a(b). 

We now introduce some notation. The sheet of b(a) with the same cut as Wa we denote 
bo (a) and the sheet joined to this via the physical cut we call 6 * ( a )- It will turn out that 
the sheets 60(a) and 60* (a) themselves are joined to other sheets via branch cuts ending at 
infinity as shown in Figure HI We denote sheets reached by traversing these cuts by bi(a) 
and &j*(a) respectively. This structure of the Riemann surface will actually be common to 
many of the functions considered in this paper and we shall use the notation for the sheets 
of such functions. 

We now proceed by demonstrating that the analytic structure of b(a) is indeed as outlined 
above and compute the asymptotic behaviour of 6(a) on each sheet using the saddle point 
equations for 6. With this knowledge we can compute both the A and B^ resolvents. 

2.1 The function b(a) below the physical cut 

The saddle point equation of B,[A] for b may be written as, 

\(W B (b) + W B .(b)) + l -(ao(b) + 00.(6)) = V'{b), (13) 

which using the fact that a (6) shares the same physical cut with Wb, (see [37]), we can 
rewrite as, 

W B (b) + ao.(b) = V'(b). (14) 

Since a and 6 are inverse functions, this may be written as, 

W B (b(a)) + a = V'(b(a)), (15) 

which using the known asymptotic form as 6 — > 00 of W B (b), W B {b) ~ 1/6, allows one to 
explicitly compute the asymptotic form of 6(a) as a goes to infinity on each sheet below 
the physical cut. In particular if V is of order p, we see that 6 has a p — 1 order branch 
point at infinity. We therefore know "half" the analytic structure of 6(a). We say "half" 
the analytic structure as we know 6(a) has a branch cut corresponding to the physical cut of 
Wa{o), however we have yet to uncover the form of the Riemann surface behind this cut. We 
therefore have the situation illustrated in Figure HI in which the other side of the physical 
cut is still unknown. 



7 




Figure 4: The structure of the Reimann surface of b(a) below the physical cut. The 
structure of the Riemann surface above the physical cut is left undetermined by the 
saddle point equations for b. 

2.2 The function b(a) above the physical cut 

To determine the analytic structure of the function on the other side of the physical cut we 
must utilise the saddle point equation for a, which for general Q is, 

(W A (a) + W M {a)) + Q(b(a) - W A (a)) = V' A {a) (16) 

which may be written, 

(2 - Q)W A (a) + b *(a) + (Q - l)b(a) = V' A {a). (17) 

Using ( fTTj) allows one to obtain the asymptotic behaviour of b(a) as a — > oo, on any sheet 
above the physical cut in terms of the asymptotic behaviour of b(a) on the sheets below the 
physical cut, where by up and down we mean the vertical direction in Figure HI Given the 
asymptotic behaviour of b in each sheet as a — > oo, an algebraic equation for b can be written 
down in which each coefficient is a symmetric function of the value of b on each sheet. More 
explicitly, suppose that b has n sheets, we make the ansatz that it arises as the solution of 
an algebraic equation of order n. Furthermore, it must have the form, 

tp n—p—2 \ p n—p—2 

J2 biAb^ + .-. + llh 6* = 0, (18) 
i=0 i=0 / i=0 i=0 

where we recall that bi is the function b(a) restricted to the ith sheet below the physical 
cut with 6 being the sheet containing the physical cut and b iJf is the ith sheet above the 
physical cut with 6 * containing the physical cut. By introducing p we have allowed for 
the possibility of there being a different number of sheets below and above the physical 
cut. The coefficients of (|18|) are symmetric combinations of the values of b on the different 
sheets. If this is to be an algebraic equation then these symmetric combinations should be 
polynomials in a. Since we know the explicit asymptotic expansion of bi and bi* we can 
compute the asymptotic expansion of each of the coefficients in (Tl~8]) . One can then see 
explicitly that either their asymptotic expansion terminates, giving only polynomial terms 
in a, or undetermined coefficients appear at positive powers of a, related to the coefficients 
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in the large b expansion of Wb{V). The requirement that negative powers of a vanish gives 
relations among these coefficients. 

This algebraic curve corresponds to the curve found via loop equation methods [3TJ and 
is known as the spectral curve. The spectral curve forms part of the initial information 
required to begin the topological recursion [29] which allows for the computation of all large 
N corrections. 

To generalise (fTT|) to allow the potentials to differ, as is required to solve (JSJ), for which 
(jHJ) is a special case, we simply redo the analysis of the saddle point equation leading to 
(j!4p . for each matrix £?W. Associated to each matrix £?W will be a function &W(a) and the 
potential for the matrix fiW w iH determine the analytic structure and asymptotic behaviour 
for &W(a) on each sheet below the physical cut. Once this information is known for each 
b^'(a), we then use the saddle-point equation for a to determine the full analytic structure 
and asymptotic behaviour of each (a) on all sheets. 

To motivate the form for the saddle point equations for a it is instructive to briefly 
consider the case of Q = 2. In this case we have three matrices, A, B^ and B^ 2 \ Having 
obtained the behaviour of the functions b^(a) and (a) below the physical cut using the 
saddle point equation for and b^ 2 \ we may find b^(a) and U 2 \a) above the physical cut 
using the saddle point equation for a which generalises (jl Tj) . 

(W A (a) + W M (a)) + (tfV) - W A (a)) + ($\a) - W A (a)) = V>), (19) 
which can be written in two ways; 

b$(a) + b$\a) = V A (a) (20) 

or 

bS?(a) + b£\a) = V A (a). (21) 

In the general case when we have matrices A and B® where i G {1, . . . , Q} all with different 
potentials one has a set of equations indexed by j, 

(2 - Q)W A (a) + bg\a) + $>«(«) = V*(a), (22) 

from which one can obtain the value of b^ on any sheet and therefore construct its algebraic 
curve. 

In the following, when working with a fixed Q and non-identical potentials, we will often 
want to avoid the bulky notation for each of the matrix degrees of freedom in (jSJ). In 
such cases the function &w(a) will be denoted by a lower case of the corresponding matrix 
in the action. For example if Q — 2 and we relabel -B^ 1 ^ — > X and B^ — > Y then the 
associated ftW(a) functions will be denoted x(a) and y(a) respectively. 

2.3 A Digression 

At this point it is useful to follow a slight digression which will produce some results useful 
in later sections of this paper. We have argued that the M resolvent for the matrix model 
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(jnj) will correspond to a free spin boundary state, however naively one might imagine we 
could make do with a simpler construction in which the free spin boundary is given by the 
A resolvent of the modelp 

Z = J [dA] [dX] [dY] exp (-NTr[V{X) + V{Y) + 3 -A 2 - {X + Y)A]j . (23) 

This expectation is incorrect since in the above model any interaction between X and Y 
is mediated by A. Hence A actually plays the role of domain walls between regions of two 
different vertex types. With this in mind, the A resolvent should actually correspond to 
placing a domain wall on the boundary. In the scaling limit we would expect this to scale 
to a fixed spin boundary state. 

Since the potentials are identical we can rewrite this model as and use the above 
equations for the case Q = 2 to check that this is indeed the case for the cubic potential 
gV(x) = x 3 /3 — x 2 — 3x. First using ffl5|) we obtain the asymptotic behaviour of b as a — > oo 
on each sheet below the physical cut, 

60(a) = ^ + l + 2^ + | + 0[a- 3 / 2 ] (24) 

Jga 2a 



61(a) = -v^+l-2^ + -^ + 0[a- 3/2 ]. (25) 



We then use the equation ( 1201) with = b^ = b to compute similar expressions for the 
asymptotic behaviour on the sheets above the physical cut. Substituting these into (I18p we 
obtain, 

6 4 - 2gab 3 - a 3 g 3 + 6 2 (-10 + g + a 2 g 2 ) + 6(-a(-10 + g)g + K ) + S + aSx = (26) 

where K,i and Si are yet to be determined constants. We contrast this with the standard 
spectral curve (j4]), which clearly differs. We find exact agreement if we shift a — > (a + b)/g. 
This shift is explained in the next section where we relate our construction to the standard 
spectral curve. However, it is interesting to pause and consider the physical interpretation 
of this shift; if we solve for a(b), clearly it only changes the result by some analytic non- 
universal constant. However, solving for b(a) and hence the A resolvent, it corresponds to an 
additive renormalisation of the boundary cosmological constant by the disc function. Letting 
the change of variable be a = (a' + b')/g and b = b', we see that b'(a') = b(a' — -b(a' + •••)) 
which one can interpret as subtracting disc amplitudes from the boundary. 6nly havine; 
done this do we obtain agreement with the known spectral curve (jl]). Algebraically we see 
it corresponds to the shift of a necessary to cancel the highest order b term. The lesson to 
learn here is that the resolvent corresponding to M may require an additive renormalisation 
of this form. 



2 Note we have chosen the coefficient of A 2 such that after integrating out A we return to the standard 
form of the action. 
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Figure 5: The full analytic structure of x(a). If the potentials V\ and V2 are of order 
p and q respectively, then the number of sheets below the physical cut is p— 1 joined 
by a p — 1 order branch cut and above the sheet we have q sheets joined by an order 
q — 1 branch cut. All branch cuts besides the physical cut extend to infinity 

2.4 The M Resolvent 

We are now in a position to begin writing down the saddle point equations for (J7j). However, 
first we must apply the method in Section 12.11 to determine the analytic structure of the 
function m(a), corresponding to the M resolvent of (J7J), below the cut. It is straight-forward 
to see that m(a) ~ a as a — > 00 and hence has only a single sheet below the physical cut. 
Furthermore, using the result of Appendix |A] we have that mo (a) = a + Wa(g>)- We can now 
write the saddle-point equations as, 

xo*(a) + y (a) = V' A (a) - a (27) 
2/ *(°) + x o{a) = V' A {a) - a (28) 
m *(a) + x (a) + y (a) = V' A {a) + W A (a), (29) 

where we have introduced x , x *, Ho an d yo* for ([7]) analogous to 6 ^, and b^} for 

The equation for m * (a) is still not particularly useful as it expresses mo* (a) in terms of 
Wx(a) which is unknown. To rectify this, we now express Wa{o) in terms of x(a). First we 
show that x(a) has the analytic structure shown in Figure [5) We already know that x(a) has 
p — 1 sheets below the physical cut, joined by an order p — 1 branch point. Above the physical 
cut we see that it differs from y(a) below the physical cut by an entire function. Hence the 
number of sheets of x(a) above the physical cut equals the number of sheets of y(a) below 
the physical cut. Furthermore, we note that the function XX=o x i( a ) possess only a single 
branch cut, identical to the physical cut of Wa(cl) and therefore it differs from Wa(o) only 
by an entire function. This is also the case for the function YHZoVii. )- Using the result of 
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Appendix |A] we find that 



p-2 q-2 , 

En— 1 - X > , v , v ''n 



£a; i (a) = W A (a)--fi and ^>i(a) = W^a) - -*=i (30) 
We therefore can write 

Vl(o) + i i^f- + tj jf) + i (|>(a) - + £ f - WbW ] • (31) 



:0 1 = 



m *(a) 



Although we could obtain the spectral curve for x(a) and ?/(a) from (|27l) and ( J28l) re- 
spectively it is easier to recall that by integrating out both the A and M matrix in (JTj) we 
obtain the standard two matrix model to which the standard analysis may be applied. We 
reproduce this analysis here, since we want to express the spectral curve in terms of x and 
a rather than the usual x and y. In particular, recalling that we may realise the (p,p + 1) 
model by setting V\ — V-i = V with V of order p, one gets 

f N 

Z = J [dX)[dY] exp(—Tc[U(X) + U{Y) - XY}), (32) 

where U{X) = gV{X) - \X 2 and g' 1 = - (1 + (7 - 1) _1 ). Recall that the spectral curve 
for the resolvent Wx(x) = ~(U'(x) — if)(x)) is E(x,ip) = where E is given in (j4]) in which 
P(x, ip) is a polynomial in x and ip with undetermined coefficient^! apart from the coefficient 
of x p ~ 2 ip p ~ 2 , where p is the degree of the potential V. A similar equation holds for Wy{y)- 
Given that a(x) = V'(x) — Wx{x) we find the spectral curve for x(a) and y(a) has the form, 

E(x, ga — x) = and E(ga — y,y) = 0. (33) 

This explains the shift found in the preceding section for the "domain wall" boundary state. 
From the above equation we can compute both x(a) and y(a) and hence, via (13T|) . compute 
m(a). Finally, also setting V\ = V2 = V with V of order p in the above saddle point 
equations, we may write, 

x *(a) + x (a) = V^(a) - a, (34) 
m (a) = I a + — — J + Xj(a), (35) 

where we have written m(a) on its 0th sheet rather than on the 0* sheet. Although ( 135]) is 
sufficient to compute the scaling limit of m(a), it is worth noting here that ( 13"5]1 can be used 
to find the spectral curve associated with the resolvent Wm by the same method one uses 
to obtain (fTj 



3 The undetermined coefficients being identicai to those appearing in (TT8|) . 
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2.5 An Example 



We give a specific example in order to illustrate the general equations introduced above. We 
again consider the case of gV{x) = x 3 /3 — x 2 — 3x. Using (13"5|) together with the asymptotic 
expansions of Xi(a) as a — > oo, which should be noted are identical to fl24l) and fl25l) . we can 
compute the asymptotic behaviour of m(a) on each sheet above and below the cut; 

m (a) = a H hO(cT 2 ), (36) 

a 

m Q *(a) = (l + ^) a -2 v ^-2-^ + 0(a- 3 / 2 ) (37) 



m u (a) = {l + g)a + 2 v ^E-2 + ^ + 0{a- 3/2 ) (3* 



5-a 

m 2 *(a) = (l + 2^)a-4-- + 0(a- 2 ) (39) 

in which the sheet m lJr .(a) and 777.2* (a) are connected by a cut of finite length. Substituting 
into a formula for m(a) equivalent to ([TBI) we obtain the spectral curve for m as a 4th order 
polynomial in a. This differs from the third order polynomial found for the spectral curve 
in [3"Ul |2"B] , however we again find we can make an additive renormalisation of the boundary 
cosmological constant m — > m + (1 + g)a such that we may remove the 4th order term in a. 
Doing this we obtain, 

q a 2 ( 1 ,\ „ „ (2 + q)m 2 2m 3 m 4 . , 

a 3 [l + m+ -m 2 ) + S ia + S + K m + ^— ^ + — + — = 0. 40 



Rewriting this in terms of the resolvent and comparing to the spectral curve in [3"Ul |2"B] we 
see that it contains the exact same terms, many with the exact same coefficient. We don't 
expect an exact agreement as we expected our construction to differ by non-universal terms. 
However, we should expect the basic analytic structure to match, and indeed in both cases 
the algebraic curve described by the above equation is a three sheeted Riemann surface with 
only finite branch cuts between sheets; which is quite different from the original Riemann 
surface for x(a). The curve fl4U|) is the spectral curve for the free boundary state and could 
in principle be used as the input for topological recursion [29]. Finally, we should stress that 
although it was necessary to perform the additive renormalisation of m to obtain agreement 
with [30j [26], this will not actually be necessary to obtain the correct scaling limit, in the 
next section. 



2.6 Scaling Limit 

In order to compare the matrix model calculation to the continuum calculation we must take 
a scaling limit of the above results. As was stated in the preceding section it is unnecessary 
to know the spectral curve for m(a); we can simply compute the scaling limit for x(a) and 
use (135]) to find m(a). 

In Daul et al. [3S], in which the method of orthogonal polynomials was used, it was shown 
that the two matrix model possess a scaling limit corresponding to the (p, q) minimal model 
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coupled to Liouville theory, when U(X) takes a particular critical form. Unfortunately, a 
formulation of our current problem in terms of orthogonal polynomials is not known. Instead 
we have saddle-point equations from which the spectral curve may be derived. 

For a given potential the scaling limit may be computed from the loop equations in the 
following way. First we fix the undetermined constants appearing in the spectral curve for 
x(a). It was argued by Eynard [31] that the condition necessary to fix these constants corre- 
sponds to the requirement that the spectral curve is of genus zero. Happily, this requirement 
means that the curve possess a rational parameterisation, i.e. there exists rational func- 
tions X : C — > C and A : C — > C such that the points on the spectral curve are given by 
(X(z),A(z)) with z e C. 

The form of the rational parameterisation, when it exists, for the curve E(x,y) = is 
known [31] and is given by, 

X{z) = - z Y J Vnz\ (41) 

71=0 

y(z) = zJ2vnZ- n , (42) 

n=0 

where we have already specialised the rational parameterisation found in [31 J to the case 
of identical potentials. Given the relations ( 133]) we then have gA(z) = X(z) + y(z). The 
dependence of the constants r\ n on the potential can be found by requiring the rational 
parameterisation reproduces the known asymptotic behaviour of x(a) as a — > oo on each 
sheet. Alternatively, explicit expressions for the constants can be found in [31]. As was 
noted in [31] , the above parameterisation first appeared implicitly in [38] via the orthogonal 
polynomial approach. 

In the orthogonal polynomial approach the matrix integral is written in terms of the 
norms of a set orthogonal polynomials p n (x). To compute the norms one needs to analyse 
the operators implementing multiplication and differentiation with respect to x. In [38] 
such operators denoted X, Xp n (x) = xp n (x), and P, Pp n (x) = d x p n (x), respectively were 
introduced and it was shown that they could be written, in the planar limit, as Laurent series 
in a single parameter z = e^Q Identifying our parameter z with the parameter z in [38] we 
see X(z) coincides exactly with X(z), i.e in our notation X(z) = X(z)\^ Furthermore, it 
was argued in [38] that P and the resolvent Wx coincide up to non-universal terms and so 
y(z) will also equal P{z) up to non-universal terms. 

For the critical potentials appearing in [38] there exists a critical value for g, denoted g c , 
such that the spectral curve possess a singular point z c where the first p derivatives of X 
and y vanish|§ The scaling limit is then defined as the blow-up of this point [2H] . 

Note that we do not require the first q derivatives of y to vanish. This is not possible for 
the following reason; the minimal potentials in [38] ensure that there exists a single point on 

4 In particular see equations (2.16) and (2.17) in [55] . 
5 See equation (2.10) in [38]. 

6 In the case of the Ising model, the critical potential corresponds to the choice, U(x) = — 3x — \x 2 + ^x 3 
and g c = 10. 
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the curve at which enough extra zeros accumulate, or equivalently enough derivative of X 
and y vanish at a branch cut to give a higher order scaling limit. However, in the symmetric 
case we have y(z) = X(l/z). Suppose z c ^ l/z c , this would mean if p derivatives of X 
vanish at z c then p derivatives of 3^ would also vanish at l/z c contradicting the preceding 
property of the minimal potential. We therefore conclude that z c = ±1 and therefore that 
y(z c + 8z) = X(z c — 5z) for small 8z, showing that y has only p derivatives which vanish at 
z c . 

We now blow up the curve to obtain the scaling limit by setting z = z c + n(g — g c ) u ( and 
g = g c + e 2 A, and letting e — > 0, 

X(z) = X{z c ) + e 2 ^P{Q + e 2 ^Q x {Q..., (43) 
y(z) = ^ c ) + e 2 ^(-l)Pp(C) + e 2 ^Qy(C) + .... (44) 

Here v = 1/ (p+q — 1) and q = p+ 1 and k is some unimportant normalisation constant. Note 
that in [38] the scaling limit of the curve (X(z), P(z)) was founc0 which by comparison to 
the above result allows us to conclude -P(C) oc T p (() and thus T p (—() = (— l) p T p ((), which 
we have used to simplify (jUJ. We may invert (jl5|l by setting ( = /i + t\i\ + 0(e 2 ) and 
solving perturbatively, to obtain, 

P(/i ) = e- 2 ^(X(z)-X(z c )), (45) 
P'M = -QxM, (46) 



which upon substitution into fj44|) gives, 

y(z) - y(z c ) = (-l) p P(/i )e 2 ^ + (Qy(n Q ) - (-1) P Q X (^)) e 2 ^ + . . . . (47) 

The above equation expresses y{z) in terms of the scaling combination e~ 2up (X(z) — X(z c )) 
which is identified with the boundary cosmo logical constant. Note that the leading order 
term in the above equation has a scaling dimension different from the known continuum 
result. This is a well known issue in the matrix model approach; the leading order terms are 
"non-scaling" and the universal physics is contained in the terms at order e 2uq . Using the 
fact that the universal terms of y{z) should coincide with those of P{z) appearing in [38] 
we conclude that Q(Q = Qy(() — (— l) p Qx(C) ^ ^<?(C)- This leads naturally to choosing 
the parameterisation_| ( = cosh(-^=). Hence for odd p we have, 

X(z) = X{z c ) + ^P{Q + ... 1 (48) 
A(z) = A(z c ) + e 2 ^Q(C) + • • • • (49) 

The case of even p will be addressed later when considering a boundary magnetic field. It 
is now useful to introduce the scaling quantities, X (a) = e~ 2up (X(z) — X(z c )) and A(a) = 
e~ 2l/q (A(z) — A(z c )), so we may write the scaled curve as, 

X(a) = A x cosh(^) A(a) = A a cosh(^), (50) 
Jpq Jpq 



7 See equations (5.8) and (5.9) of [38] . 

8 The coefficient of a is chosen to agree with the literature. 
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where A x and A a are constants. To find the scaled version of the curve corresponding to the 
Wm resolvent we must find an expression for the scaled version of Xi(a). Given a single value 
for a there exist different values of x(a), Xj(a), corresponding to the different sheets, hence 
in the scaling limit we must find the values of a that give the same value of a while changing 
the value of x. Since we have A(a) oc cosh(-^=), we see that a — > a + 2my^|, where fiGZ, 
moves between the sheets of x(a). We therefore have from (|35|) . the scaled spectral curve, 

M(a) = 4*Xjcosh^-^= ( a + 2ni yf)) ' ^( a )= A - cosh (^7=)- (51) 



pq 

Note in particular that the additive renormalisation of m introduced in the previous section 
would not affect this result, as it would only contribute sub- leading order terms. At this 
point we also want to draw attention to the similarity between this relation and the Seiberg- 
Shih type relation appearing in the introduction. In fact we can make this similarity more 

manifest by making the change of variable a — > a — i./-(p — 2), giving, 




M(a) = A x >J cosh ( ( cr + nij~ ) ) , A{a) = A a (-l) q+1 cosh(^). (52) 



By using the expression for the Chebyshev polynomial of the second kind, Uk(cos(x)) 



Yln=-k 2 ex P (i nx )i we see that, 



M(a) = A,(-l)« +1 ^_ 2 (cos(-))cosh(^), A(a) = A (-l)« +1 cosh(^), (53) 

q Jpq Jpq 



where we have used the fact that p + 1 = q. To make contact with the usual form of the 
continuum disc amplitude it is necessary to rescale the matrix M appearing in the resolvent 
by (— l) 9+1 £7p_2(cos(-)). This rescaling was already necessary in the work [301 EH]- Once 
rescaled, the scaled spectral curve is exactly the expression one would expect for the disc 
with a (l,p — 1) boundary state. We therefore identify this state as the free spin state of 
the (p,p+ 1) model. 

The expression for Ai(a) in (I52p is very reminiscent of the Seiberg-Shih relations. How- 
ever, there are some key differences. Firstly, this relation applies to resolvents rather than 
the non-marked disc function appearing in the original relation. Secondly, this relation arises 
as the scaling limit of a discrete version of the relation fl35|) . Unfortunately we have only 
constructed the (1,2?— 1) boundary state; this is a state that has been constructed before in 
the matrix model [2H] . Our construction is nonetheless interesting, as it provides an explicit 
spectral curve from which this resolvent arises and therefore it is much more obvious how one 
might begin to apply the topological recursion [29] to our construction. Indeed it appears 
simple to generalise the discrete expression (1331) so that it produces a scaling limit consistent 
with other boundary states, 



m {a) 



a + ^i +5> 4 (a), (54) 



p 
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where 2 < I < p — 1. 

This construction sheds light on the results in [25l [26] in which it was observed that the 
Seiberg-Shih relation no longer appears to hold for worldsheets of higher genus. It would now 
seem there is truth to the Seiberg-Shih relation; it gives the relation between the spectral 
curves of the different boundary states. However, we conjecture that to correctly compute 
the higher genus amplitudes one should use the appropriate objects defined on these algebraic 
curves in [29] rather than naively applying the Seiberg-Shih relation to every amplitude. 



3 Boundary Magnetic Field 

We now consider the situation where "spin up" and "spin down" fields will be given different 
weights when appearing on the boundary; this should correspond to a boundary magnetic 
field. Particularly, we aim to generalise the boundary magnetic field studied for the Ising 
model in [5UJ [25] for all (p,p + 1) models. Consider the action, 

Z = j [dMpApxpr]e- WTr ^ (55) 

where the potential Vh has the form Vh(x) = ^-x 2 + Y^T=i ^t xk - Integrating out M and A 
results in the two matrix model with critical potentials independent of h. Hence, ones sees 
that there is no h dependence in the bulk. However, the M resolvent will be h dependent. 
From a similar argument as presented around ([6]) and Figure [3], one then expects that the 
M resolvent will flow to the same boundary conditions as the e h X + e~ h Y resolvent. This 
resolvent is usually referred to as the free spin disc function with a boundary magnetic field 
h [30]. Indeed, one sees that the "spin up" fields (corresponding to X's) adjacent to the 
external M k vertex in Figure [3] are weighted with a factor e h while the "spin down" fields 
(corresponding to Y's) adjacent to the external M k vertex are weighted with a factor e~ h . 

We may introduce new variables R = e h X and S = e~ h Y in order to return the action 
in (|55|) to the form in (J7J), 

Z = J [dM][dA][dR][dS]e~ NT ^^ +V -^^^ (56) 

where Vh(x) = Vh(e~ h x). Clearly the saddle point equation for m(a) remains the same, while 
the equations for r(a) and s(a) will differ from x(a) and y(a). Since we will express m(a) 
in terms of r(a), we must find a rational parameterisation for the curve r(a). This can be 
accomplished by the fact it is possible to relate r to x and s to y. Upon integrating out M 
and A to obtain the standard two hermitian matrix model, 

f N 

Z = J[dR][dS]exp(—Tr[U h (R) + U^ h (S)-RS]), (57) 

with Uh(x) = gVh(x) — \x 2 and g^ 1 — - (1 + (7 — 1) _1 ), we immediately know the form of 
the spectral curve |3T] . 

E h (r, s) = (U' h (r) - s) {U'_ h {s) - r) + gP h (r, s) = 0, (58) 
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where s(r) is related to the resolvent by Wr{t) = ^{U' h (r) — s(r)). Noting that U' h (x) = 
e U'{e x), we can write, 

E h (r, s) oc (U'(e- h r) - e h s) (U'(e h s) - e~ h r) + gP(e- h r, e h s) = 0, (59) 

where we have rewritten by extracting the dependence on the magnetic field from the 
arbitrary constants. Finally given a is related to s in the above equations via V^(r) — a{r) = 
-{JJ' h {r) ~ s ( r ))> we can wr he, 



E (e~ h r, ge h 



-U' h {r) - K(r) + a(r) 



0, (60) 



where E is the spectral curve for the two hermitian matrix model in which both potentials 
are U. Since the curve E in the above is independent of h, we can relate the rational 
parameter isation at h ^ to the h = case, 

H(z;h) = e h X{z) (61) 
A(z;h) = Vl{e h X{z))-e- h [V^X{z))-A{z)] 
= e- h W(X(z)) - V^X{z))) + e- h A(z). 

The saddle point equations for m can now be written as, 

/ t \ *~ 2 

M{z; h) = A{z; h) + e h ^ + e h ^ X{U\ 0), (62) 

\ p ' i=0 

where corresponds to a point on the zth sheet of x(a), i.e. given z we can find Zi with the 
same value of a but differing x. This represents a generalisation of the Seiberg-Shih relations 
to the case of non-zero boundary magnetic field. 

We can now also address the issue of the scaling limit in the case of even p. Consider 
flfJTj) when h = iir/2, we have that, 

A(z; itt/2) = -~{gA{z) - 2X{z)) = --(y(z) - X{z)). (63) 
9 9 



If we substitute in the scaling form for X and y, i.e. f H3|) and OH]), we find that to leading 
order, 

A{z;i*/2) = —f^QiQ (64) 
9 

and hence the scaling limit of Ai(z; in/2) equals the sum over X{ti\ 0) up to a normalisation. 
Hence for even p and h = in/ 2 we have that the scaling forms of A(z; in/ 2) and M.{z\ in/ 2) 
are equal to the scaling forms of A(z) and M.(z) for odd p with h — 0, the analysis of Section 
I2.6l then applies directly. We therefore conclude that for even p the free spin boundary state 
is created by the X — Y resolvent. 
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3.1 Scaling Limit 



We now consider the scaling limit for general values of h for odd p. We expect even p to 
be identical, however it requires a slight modification of the algebra. Recall the expression 
flfJTj) . Given that we know X(z) and A(z) scale as fHHl) and (HDl . we can obtain the scaling 
form for A(z; h), 

A(z; h) = a c + e~ h (Vax c ) - V^x c )) + e~ h (V^'(x c ) - V£{x c ))A x CO sh(^)e 2 ^ + 



+ e~ h ((V£(x c ) - V>\x c ))F{o) + A a cosh^)) e 2 ^ +1 > + . . . , (65) 



where F is an unknown function defined implicitly by X(z) = x c + A x cosh(-^|=)e 2 ^ p + 

F(a)e 2u(j}+1 ^ + 0(e 2u ( p+1 }). The crucial point to notice about this scaling relation is that for 
h 7^ 0, a — (a c + e~ h (V^(x c ) — Vq(x c ))) = 0(e 2up ), whereas for h = this term vanishes 
and the leading order behaviour of a is controlled by the next-to-leading order term. This 
property affects the sum over the sheets Xi(a) in the expression for m(a); the value of x 
on its different sheets, Xi, are obtained from the values of a which map to the same a, but 
differing values of x. By considering the expression for a in ( )65l) . we can see given a point 
on the curve a, that the points on the curve a' n given by, 



' a + 2in^+e 2v K n {a) + 0{e Au ), if h ^ 



< = { V r (66) 

a + 2in^+0{e 2l/ ), if h = 0, 



map to the same values of a to order e 2u( - p+1 \ where K n is given by, 

A 'Vr* ] J ) = ( ' " ( + ° KM -KM ' (6?) 

where cr n = a + 2in./-. Substituting the h ^ case into m we obtain, 



M(z-h) = m c + B(h)A x cosh(^)e 2 ^ + Q(h) 



cosh(^) 



e 2^(p+l) 



^ cosh (^), (68) 



where Q(h) = e~ h (V/ l '(x c ) — Vq(x c )) + (p — l)e h . Eliminating a from fl68l and fl65|) we 
find that the first non-universal term exactly reproduces the disc amplitude for the identity 
boundary condition. This extends the picture found in [30l 126] for the Ising model; for any 
(p,p + 1) model we find the boundary condition corresponding to a magnetic field has two 
distinct fixed points; one corresponding to free spin for h = and the other to fixed spin. 

Finally, we can also comment on the dual brane found for certain values of the magnetic 
field in [26]. It is clear from the scaling form for m that for values of h satisfying Q(h) = 0, m 
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must scale with a non-physical dimension. For the case p = 3 this is precisely the dimension 
associated with the dual brane, as was necessary in (26] to obtain a non-trivial limit. Finally, 
substituting the expression for m into a in the case Q(h) — results in precisely the dual 
brane amplitude. Although this qualitatively reproduces the results in [26], the precise 
value of h for which the dual scaling limit occurs in j2B] was h = m/2 which differs from 
that found here. A possible explanation of this can be found in our earlier discussion of 
the renormalisation of the boundary cosmological constant by the disc function, which was 
necessary in order for the spectral curve for a(m) to match with that found in [26]. This 
renormalisation took the form m — > m + rja, which modifies (|68|) by replacing, 



Requiring that we reproduce the result of [26] we find 77 = g c (p — l)/2 + 1. For p = 3 this 
reproduces exactly the additive renormalisation used in Section 12.51 

4 Discussion 

In this paper we have found a number of generalisation of previously known results. Firstly, 
we have shown that by simply introducing an extra matrix with a Gaussian potential into 
the two-hermitian matrix model we are able, by computing the resolvent of this new matrix, 
to variety of new boundary states corresponding to a non-zero boundary magnetic 

field. We were able to show that for the (p,p + 1) minimal string, the (l,p — 1) state 
corresponds to the free spin boundary and that when a non-zero boundary magnetic field is 
introduced all discrete states flow to the fixed spin boundary state in the continuum limit. 
This generalises previous work in [301 ES] • We found also the surprising result that for even 
p the free spin boundary is created not by a X + y-like resolvent but by a X — Y resolvent. 

In order to compute the resolvent corresponding to the new state we have generalised the 
methods of [36J to the case of non-identical potentials and shown how the structure of these 
equations leads to Seiberg-Shih like relations between the various boundary states. The key 
difference between the known Seiberg-Shih relations and the ones appearing in this paper 
is that ours are valid away from the continuum limit. Indeed, it is trivial to construct the 
spectral curve associated to the new boundary states. This leads to an interesting possibility 
for resolving the tension between the success of the Seiberg-Shih relations in the case of 
disc functions and their apparent failure for more complicated amplitudes; the Seiberg-Shih 
relations relate the algebraic curves associated with the different boundary states. If this 
is true it would mean that the identity brane does indeed contain all the information for 
constructing the other branes, however the higher genus amplitudes should be constructed 
by applying the topological recursion relations [20] rather than applying the Seiberg-Shih 
relations directly. One obstruction to this program is the necessity to better understand the 
additive renormalisation we found necessary to bring the spectral curves found here in to a 
form found previously. Although this renormalisation was not necessary to obtain agreement 




(69) 
(70) 



(p - 1) cosh(/i) + —(2(1 - 77) + g c {p - 1)) sinh(/i). 
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in the scaling limit, it is unclear whether the renormalised or non-renormalised curve should 
be used for the topological recursion. We hope to pursue this possibility in future work. 
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A Appendix 

We prove the following useful result: 

Lemma 1. Given a potential V(b) = X^fe=i %^ fe with n > 1, then the asymptotic behaviour 
ofb(a) as a — > oo on a sheet below the physical cut, as determined by the equation Ws(6(a)) + 
<j>a = V'(b), is 6(a) = (0/t n+1 ) 1 /" a ^ - ^ + 0{a~ l / n ). 

Proof. Recalling that the asymptotic behaviour we want is determined by the equation 
PUe(6(a)) + 4>a = V'(b), we may write, 

n 

Y,tk+ib k = 0a + O(a- 1 ), (71) 

fc=0 

where we have used the known asymptotic behaviour of Wb(P) ~ 1/6 on its physical sheet. 
Making an ansatz for the asymptotic form of 6(a) of, 

— oo 

6(a) = a m a^ (72) 

m=l 

we have our desired result by substituting this into ( 17 ip and comparing coefficients of a. □ 
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